Non-analyticity in scale in the planar limit of QCD 
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Using methods of numerical Lattice Gauge Theory we show that in the limit of a large number 
of colors, properly regularized Wilson loops have an eigenvalue distribution which changes non- 
analytically as the overall size of the loop is increased. This establishes a large- TV phase transition 
in continuum planar gauge theory, a fact whose precise implications remain to be worked out. 



Intuitively, one expects parallel transport round a 
closed curve in four-dimensional SU(iV) pure gauge the- 
ory with 9c p — to be close to identity for small 
curves and far from identity for large curves. In this 
letter we make this idea concrete and find that small and 
large loops are separated by a large- N phase transition. 
The possibility of a new type of non-analyticity entering 
SU(AT) gauge theory in the 't Hooft limit p] N — > oo has 
preoccupied researchers for a long time; here we present 
a class of examples where this phenomenon occurs. The 
numerical evidence is sufficiently convincing to view the 
effect as an exact property of the large- TV limit. We start 
presenting a simple lattice result and let it lead us to the 
above conclusion. After that, we put the result in his- 
torical context, mentioning some of the main papers this 
work is related to. 
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FIG. 1. Histogram of eigenvalue angles for a Wilson loop. 



The gray line in Fig.[T]is the histogram for the angles of 
the eigenvalues of a square loop in SU(29) lattice gauge 
theory with side length equivalent to about 0.54 fermi 
in QCD units and viewed at a resolution (thickness) of 
about 0.15 fermi. The black line is the eigenvalue-angle 
density determined by the heat-kernel function (HK) for 
SU(29) which depends on a single parameter t, to be 
defined later (in the plot, t w 4). The two curves cannot 
be distinguished in the figure, so the HK approximates 
the data well, after adjusting t. The data is sufficiently 
accurate to check whether the HK might provide an exact 
description. This is definitely ruled out by a \ 2 analysis. 

Some definitions are in order now. The Wilson loop 



matrix associated with a closed spacetime curve C is 
W r (C,x,s) = 7>exp (i^A r fl (y,s)dy l A € SU(iV), 

where r denotes an SU(iV) representation, a; is a point 
on C, and s > denotes a "smearing parameter" of di- 
mension length squared (y/s determines the thickness of 
the loop). Smearing is required to make all W r (C,x,s) 
finite SU(iV) matrices with operator-valued entries. The 
gray line on the previous plot shows 
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The 0i(C,s) are angles locating the eigenvalues of 
W/(C,x,s) (/ denotes the fundamental representation) 
on the unit circle; they do not depend on the choice of 
x. After averaging, also the dependence on the location 
and the orientation of the loop drops out. With the CP 
violating 9 c p parameter set to 0, pn is invariant under 
9 -> 2tt - 6. 

Smearing is defined as follows: One starts with five- 
dimensional gauge fields on ]R 4 x M + ; the smearing pa- 
rameter s lives on the M + . The usual quantum fields are 
denoted by B^(x) and reside on the R 4 -boundary. The 
Af(x, s) are defined for s > by 

F f tS = Df^Fl v with A{ l (x,s = 0)=Bf l (x). 

The 5D gauge freedom is reduced to a 4D one by 
A((x,s) = 0. At s > 0, all divergences coming from 
coinciding spacetime points in products of renormalized 
elementary fields are eliminated by a limitation on the 
resolution of the observer, parametrized by s. Renormal- 
ization of the boundary quantum-fields Bj^(x) proceeds 
as usual. The definition of smearing easily extends to any 
finite UV cutoff including the lattice: replace DyF^ ^ 
by the variation of the action with the UV cutoff in 
place. Smearing extends formally to loop space, with 
C parametrized by cr, 

d s Tr(W f (C,s)) = f f^iC,*)) ^ LTl{Wf{CjS)) . 

L is the Levy Loop Laplacian appearing in the Makeenko- 
Migdal (MM) equations. Were a string representation 
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of L found, diffusion in loop space would become well 
denned and field theory and string theory could refer 
then to the same non-singular object. 

The HK (heat kernel) probability density (w.r.t. the 
Haar measure) for an SU(TV) matrix W is 
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implying {Xr(W)) — d r e"^ C2 ^ r \ The parameter t is a 
"diffusion time" and d r , C2(r) are the dimension and the 
quadratic Casimir of the irreducible representation r (in 
Fig. [I] t — 3.881). The heat kernel represents a mul- 
tiplicative random walk on the SU(TV) group manifold 
emanating from the identity. The HK single eigenvalue 
distribution 27rp]Q- K (0, t) is given by 
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The HK represents the data very well, but is not ex- 
act. One could explain this by postulating "Casimir 
dominance": Tr (W r (C,s)) d r e -CaW©(C, a ) ) with a 
r-indcpendent &(C,s). This approximation must break 
down for very large loops, where screening effects come 
in and the loop is dominated by the area term. Then, 
only the TV-ality of r should matter. However, both in 
perturbation theory and at intermediate scales (up to 2 
fermi) , Casimir dominance is known to be a good approx- 
imation. All our data is in the range (0,1) fermi because 
the large-TV transition is roughly in the middle of this 
segment. 

We proceed to make the case that this indication of 
approximate Casimir dominance can be replaced by an 
exact statement in the large- TV limit, namely that the HK 
formula is exact in a precise "large- TV universality" sense 
to be explained below. The main point is that the shape 
of the eigenvalue-angle distribution of smeared Wilson 
matrices associated with uncomplicated loops is governed 
by two opposite tendencies: random-matrix eigenvalue 
repulsion and asymptotic-freedom attraction to unity. 
This produces TV alternating peaks and valleys with a 
swing of order between them. So, most of the structure 
one sees in Fig. [I] is determined by a generic mechanism. 
Rather than the oscillations, the truly interesting feature 
is the deep well around n. For a smaller loop, this well 
would be deeper and wider. Taking TV to infinity elim- 
inates the more obvious details and leaves the essential 
features: for a small loop the valley around it flattens out 
at TV = oo and the eigenvalue density is zero there. For a 
large loop, the eigenvalue density is non zero around the 
entire unit circle. In the HK case, the role of overall loop 
size is played by t and the transition between a gap-less 



and a gapped spectrum occurs at t = 4. Large- TV univer- 
sality will be a statement about the large- TV behavior of 
the gauge theory data in the neighborhood of a loop size 
that would be critical at N = oo. 

To concentrate on the region of interest we need to 
choose an observable that is particularly sensitive to 
eigenvalues close to -1: 



N (v,C,a) = (det (e*+e-*W>(C,*)) 
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This observable generates all Wilson loop expectation 
values in the fc-antisymmetric irred. representations of 
SU(TV) and is given by £f =0 e (^" fc > ( x i sym (W f (C, s))). 

At the transition point, logOjv will have a non- 
analyticity at y = when TV = oo. Therefore, we con- 
sider the expansion 

N (y 7 C, s) = a Q (C, s) + a x {C, s)y 2 + a 2 (C, s)y i + O (y 6 ) 

and use ujn(C, s) = "'^p' ) as our new observable. 
In the HK case, with r = t(l + l/N), ^ K (y,r) 
v (jMog(0^ K (?/,-r)) obeys Burgers' equation: 
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N enters only via the viscosity, which is We may 
replace the term "large- TV universality" by "Burgers uni- 
versality". Our observable in the HK case is plotted in 
Fig. [2] for different values of TV. 




FIG. 2. Development of a jump singularity in the HK case. 



Figure [2] shows that w]J K (r) becomes a 9- function at 
infinite TV but the singular behavior develops slowly with 
TV. Burgers' equation implies that the jump in the 6- 
function causes d y (f)^(y = 0,r), which is finite and 
negative for r < 4, to blow up when t reaches r = 4. 
We expect a similar behavior in the four-dimensional 
gauge theory. To be specific, Burgers universality means, 
e.g., that the following predictions about the large-TV 
limit hold exactly in the four-dimensional gauge the- 
ory, where criticality sets in at an overall loop size 



3 



/ = L: lim 



N- 



!1± 



l=l a 



3 J_ 
2 if' 



(here, K 



— r 2 f 1 ) 

4w I 4/ 



"I ^2- 



1 -\ ^K- 

24 V 2 ' 



1.046); limjv^ooiV- 

lirri7v->.oo Wjv|; = ; = g^ 2 - Here, the variation w.r.t. the 
overall size I is taken at constant loop shape. The roots 
of On{u) are all on the imaginary axis (as a consequence 
of the Lee- Yang theorem) . Another universal property is 
that in the critical regime (around y — 0, I — l c ) these 
roots scale like N~ s . 

Our objective was to establish numerically that the 
transition occurs and that the universal predictions hold 
in the continuum limit of a lattice gauge theory with 
standard single plaquette Wilson action. The lattice cou- 
pling, traditionally denoted by /3, is determined by the in- 
verse 't Hooft coupling b = nW- As N varies, simulations 
are carried out in varying ranges of 6, all contained in 
the segment [0.348,0.380] with upper limits determined 
by the lattice volumes, V . We have carried out simula- 
tions at N = 11, 19, 29 on hypercubic lattice volumes, V, 
(12 4 ,14 4 ,18 4 ), (10 4 ,12 4 ,13 4 ,14 4 ), (8 4 ,10 4 ,12 4 ), respec- 
tively. The volume size determined the maximal b al- 
lowed in order to stay in the confined phase at infinite N. 
For finite and fixed N, at b values close to the maximal 
allowed one, sizable finite- volume effects became evident; 
they were consistent with an exponential dependence on 
the linear size of the system. The data at different vol- 
umes was thus used to eliminate results contaminated by 
finite-size effects. 

The measured observables were extracted from square 
L x L Wilson loops, with L = 1, 2, . . ., 9. Smearing was 
implemented on the lattice with a parameter S = L 2 /55 
(our convention is that capital letter symbols correspond 
to lower case symbols in the continuum). For a given 
Wilson loop, we extracted the coefficients do, 1,2 for each 
orientation and location of the loop. After averaging we 
obtained one set of coefficients for each gauge configura- 
tion. We generated 160 independent gauge configurations 
at each set of parameters (b, N, V). Consecutive b- values 
were spaced by increments Ab = 0.001 and were sepa- 
rated by 1000 passes, half of heat-bath type and half of 
over-relaxation type. Looking at autocorrelations, we de- 
termined that they typically drop by e for every 250 such 
passes. Averaging over configurations produced coeffi- 
cients which determined our estimate for u>pf(b,L), and 
their errors. The data were sufficiently closely spaced in 
b that we could use spline interpolation to represent the 
result as a continuous function of b. 

From the HK case we learned that prohibitively large 
values of N would be needed to display directly the singu- 
lar large-A behavior. We adopted instead the following 
strategy: From each ujn (6, L) we constructed a number 
tat(&, L) by solving u>^(b,L) = w™(tjv(&, L)). The re- 
quired inversion is unique. We then showed numerically 
that the convergence of tat(6, L) to ^(6, L) is rapid, like 
that of the string tension and the deconfinement tem- 



perature. Hence, LJoa(b,L) = liniAr^co u;^ K (Tjv(&, L)) = 
liniAr^oo w^ K (r 00 (6, L)) and to some universal subleading 
order in w/v(&, L) sa w™(Too(&, L)). This relation can 
be taken over to the continuum. Extrapolating to the 
continuum limit, the two variables b and L get replaced 
by a single length variable, I, the side of the loop in phys- 
ical units. We took measurements in a region in which 
the emerging Too(Z) extends on both sides of r M = 4. We 
numerically determined that ^(1) is smooth in I and in- 
vertible at I = l c where t^I^) — 4. We conclude that 
the continuum ujn(I) will develop the same singularity in 
the vicinity of I — l c as ai^ K (Too(/)) would, establishing 
the transition and its universality. 







0.4 








0.3 












-1.0 


-0.5 


0:2 
0.1 
0.0 





2 log(A r /) 



FIG. 3. Rapid convergence to N = 00 for rjv(Z) in continuum. 



The continuum functions — Vn are shown in Fig. 3 for 
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N = 11, 19, 29 by a dashed line, a solid line, and a aot- 
ted line, respectively, and rapid convergence to N = 00 
is evident. The horizontal axis is labeled by 21og(Ayi). 
At is the infinite-A critical deconfinement temperature, 
roughly equivalent to 264 MeV in QCD units. The shades 
indicate the accumulated errors, for those regimes where 
we had enough data to reliably estimate them. By that 
we mean that we had data points at least at three dis- 
tinct pairs (b, L) all corresponding to the same physical 
I. The relation between I and (b,L) was set by a stan- 
dard one-loop tadpole-improved formula, known to work 
well from other simulations. A reliable estimate of the 
continuum limit could be obtained by extrapolating in 
the lattice spacing squared and observing a linear behav- 
ior. At N = 29, only two pairs were available for some 
values of / and a continuum number was obtained by 
postulating a linear behavior, but without an error. The 
approximate linearity of the continuum functions is 
consistent with asymptotic freedom when tn(1) is viewed 
as an effective running coupling constant. The slope of 
the lines in Fig. [3] is about 0.22, while the expected co- 
efficient is about 0.29. The discrepancy has to do with 
a factor associated with smearing, which has been so far 

only calculated at tree level where the theory is conformal 

.2 

with a dependence only on the ratio y. 

Because of our indirect way to establish criticality, 
some of the large-A universality predictions we listed 
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earlier become tautological. However, there remain two 
extra checks that are meaningful: one is to determine the 
exponent 3/2 from the ratio do/ai and the other is the 
exponent 3/4 from the N~ 3 / 4 level density in the critical 
region. Logarithmic fits produce estimates within 1% of 
the expected values. This concludes our account of the 
numerical evidence for the non-analyticity. 

Having this one example of a large- N phase transition 
opens up some new questions: Is the transition physi- 
cal in the particle physics sense, that is, would one ac- 
tually see new singularities in a large-TV, narrow-width 
approximation to the S'-matrix of the theory? If the an- 
swer is positive, one might speculate that at N — oo 
one does have exactly linear Regge trajectories for high 
spin states, up to a point, beyond which, presumably, a 
perturbative behavior ensues. If the answer is negative, 
one would need to understand how exactly the S-matrix 
gets shielded from the above large- N phase transition we 
have established. In this context it is important to note 
that the non-analyticity we found does not occur in the 
Wilson loop expectation values themselves, so long as the 
number of boxes in the Young tableau of the representa- 
tion is kept finite and fixed as N — > oo. 

Many avenues for further investigations are likely to 
open up, providing some encouragement to those pur- 
suing the long quest of conquering the large-iV limit of 
QCD, albeit, perhaps, in only a semi-analytic way. 

We now turn to a review of the history of the sub- 
ject, focusing on the main contributions. We avoided do- 
ing this earlier in order to keep the presentation stream- 
lined. The large- TV phase transition we discuss was dis- 
covered by Durhuus and Olesen [5] in the context of two- 
dimensional SU(iV) gauge theory, where Casimir domi- 
nance is exact. They obtained the phase transition work- 
ing directly at infinite N, where they identified the invis- 
cid Burgers' equation as playing a central role. Therefore, 
it would be fair to refer to this transition as the DO tran- 
sition. Blaizot and Nowak [3] argued that the large- N 
universality class was controlled by Burgers' equation [3] , 
including its viscous term. This was shown to hold ex- 
actly in the HK case [5] . The generic nature of the DO 
transition, in the sense that it occurs also in three- and 
four-dimensional SU(iV) gauge theory, was conjectured 
about five years ago [6[. A test in three dimensions was 
shown to be consistent with the conjecture, but there 
was no overwhelming numerical evidence [7] . Continuum 
smearing was introduced at the same time, as a techni- 
cal device used to define "eigenvalues" of Wilson loops in 
renormalized QCD. The approximate validity of a "dif- 
fusive" viewpoint (that is, using the HK as a model) for 
the behavior of Wilson loops has been pointed out al- 
ready in 2005 [8], for the case of the gauge group SU(2). 
The exact formula for the single eigenvalue- angle density 
in the HK case at any N was derived in [3]. Casimir 
dominance has been discussed at the perturbative level 
in [10] and at the nonperturbative one it was reviewed 



by Greensite in [TTJ . The limitation on the size of a finite 
box at infinite N is studied in [12]. 

The four-dimensional test in this paper became pos- 
sible after a hardware upgrade of a computer cluster at 
Rutgers. Substantial resources were invested to produce 
a convincing and detailed case in four dimensions, as the 
question of whether a large- N phase transition can oc- 
cur in four-dimensional continuum SU(iV) gauge theory, 
as the scale is varied, has been considered repeatedly in 
the past, but without definitive evidence one way or an- 
other [T3"] . 

The relation between the large- N expansion and the 
program of dual topological unitarization [T3] was ex- 
plored a long time ago by Veneziano [TS] . The associated 
question of Regge trajectory linearity has been addressed 
by McGuigan and Thorn [16]. A more direct connec- 
tion between the infinite- N limit and zero string-coupling 
string theory has been an object of study for a long time, 
and perhaps the best known line of attack is the sugges- 
tion to start from the Makeenko-Migdal [T7] equations 
and guess a solution defined by a string theory. 
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